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1. Introduction 

Since tlie discovery of high temperature superconductivity in cuprates [Q], a tremendous 
effort has been made to uncover the mystery of this phenomenon. It is generally 
believed that the strongly correlated electron systems behaving as non-Fermi liquids 
are closely related to superconducting materials. This has caused an intense study in 
strongly correlated electron systems |], H, ||, 0]. These systems possess various 
physical characteristics which are decisively dominated by the competing interactions; 
e.g. the Coulomb interaction in the Hubbard model, spin fluctuations through the 
antiferromagnetic coupling for the super-symmetric t-J model and current-density 
correlated interaction inducing hole pairs of Cooper type superconductors in the one- 
dimensional (ID) Bariev model. The ID Hubbard model as a prototype amongst the 
strongly correlated electron systems has attracted a substantial deal of interest in the 
study of integrable quantum field theory, mathematical physics and condensed matter 
physics since it's exact solution was achieved by Lieb and Wu in 1968. Towards a 
complete understanding of the mathematical structure of the ID Hubbard model in the 
framework of the quantum inverse scattering method (QISM), a fundamental advance 
was achieved by Shastry |^ in demonstrating the integrability of the model. Specifically, 
it was shown that a two-dimensional statistical covering model of two coupled symmetric 
six vertex models provides a one parameter family of transfer matrices commuting with 
the Hamiltonian of the ID Hubbard model. The algebraic formulation with respect 



to the integrability leads to the quantum i?-matrix P, 10, 11, O] which facilitates not 



only the algebraic Bethe ansatz solution |TB[, but also the construction of the boost 



operator ||T4[ for the model. Remarkably, the Hamiltonian of the Hubbard model was 
proved to exhibit the 5*0 (4) symmetry by Yang and Zhang [jl5| (see also ||16|). Besides 
the spin SU{2) algebra, the SO {A) algebra contains the r^-pair SU{2) algebra with the 
raising operator creating a on-site pair of electrons with opposite spins. This can be 
interpreted as a localized Cooper pair. A complete set of eigenstates of the Hamiltonian 



can be obtained by exploiting the SO {A) symmetry |jT7 . 

The ID Hubbard Hamiltonian with more competing interactions may also be 
considered. Along this line, many extended Hubbard models have been constructed 
in the literature, such as a u{2\2) extended Hubbard model 0, supersymmetric 



Uq{osp(2\2)) electronic systems [|T^ and SU{N) Hubbard models |jT9[. In this paper, 
we present an alternative ID Hubbard model such that the Hamiltonian has off-site 
Coulomb interaction instead of the on-site one of the standard Hubbard model. The 
integrability of this model is verified by presenting the associated quantum i?-matrix 
which fulfills the Yang-Baxter equation (YBE). We show that the model exhibits the 
5*0(4) symmetry with new representations of the r^-pairing SU{2) algebra and the (- 
pairing spin SU{2) algebra. Moreover, the algebraic Bethe ansatz is formulated by 
means of the QISM. Though the model exhibits the same spectrum as the standard 
Hubbard model on a periodic lattice, the new quantum i?-matrix, the hidden nesting 
structure associated with an asymmetric isotropic six-vertex model and the Bethe 
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eigenvectors do distinguish this model from the standard one |T^. The essential 
differences between the two models manifest in the open lattice versions, which we will 
discuss in more depth in the conclusion. 

The paper is organized as follows. In section ^ we introduce a Lax operator 
associated with the new Hubbard model and construct a nontrivial higher conserved 
quantity commuting with the Hamiltonian. In section ^, we present the i?-matrix 
associated with the model by solving the Yang-Baxter relation. The 5*0(4) symmetry is 
verified too. In section we formulate the algebraic Bethe ansatz solutions for the model 
with periodic boundary conditions. The eigenvectors and eigenvalues of the Hamiltonian 
are presented explicitly. Section || is devoted to a discussion and conclusion. 



2. The model 



Let us begin by introducing a variant of the ID Hubbard model with the Hamiltonian 

^ = E { + + - ^)} + 7 E ^1^+1- (2-1) 

j=i i=i 
Above aj and tj are the two commuting species of Pauli matrices acting on site j, f/ is a 
Coulomb coupling constant. Above and throughout, periodic boundary conditions are 
imposed on all summations evaluated over the lattice length L. The difference from the 
standard Hubbard model is that the model ( ^.1[ ) exhibits the off-site Coulomb interaction 
instead of the on-site one. We shall see that it not only breaks the spin reflection 
symmetry but also specifies a new representation of r^-pairing SU{2) algebra and spin 
SU{2) algebra in order to maintain the 5*0(4) symmetry. To verify the integrability 
of the model, we, at first, identify a relation between the Hamiltonian (|2.1|) and the 
transfer matrix which is defined by 



t{u) = TroT{u) 



[2.2) 



with 

T{u)=LoLiu)---Loi{u). (2.3) 
The local Lax operators associated with model ( |2.1| ) have to be alternatively chosen as 



where 



qKu) ^ 



-h{u) 



-h{u) 



a-Q] 



e-Hu)p-Q+ 



-h{u) 



^-h(u)p-^- 



:2.4) 



:2.5) 



Integrable variant of the one- dimensional Hubbard model 



4 



with a parameterization ■y{u) = Wi{u) —Wsi^u) = sin(-u); a{u) = w^^u) +W3{u) = cos(m). 
We would hke to mention that the Lax operators 

L^j{u) = w^iu) + W3(M)(T|a^ + (T+aQ- + a+aj (2.6) 

L^^iu) = W^{U) + W^{u)t^T^ + T+Tq- + T+TJ (2.7) 

/o = cosh^ + aoXsinh^, (2.8) 

have been chosen the same as that for the Hubbard model [1^, |lT|. It follows that 
the Hamiltonian ( |2.1| ) is related to the transfer matrix matrix ( |2.2| ) in the following way 

In t{u) = In r(0) + H u + ^ J u'^ + ■ ■ ■ (2.9) 

above the Hamilonian H = X]j=i with the Hamiltonian density 

= Loo+i)(0)L;^.(0)Lo/(0))Lo(;+,)(0), (2.10) 
and the second higher conserved current can be given as 

L 

J = ^ Jj(j+l)(j+2) (2.11) 



with 



J3U+W+2) = - Hj^j+i) - [Hj{j+i), i^'(j+i)(j+2)] , (2.12) 

= Loo+i)(0)L;.(0)L-1(0)L-(;.^^)(0). (2.13) 

Here the prime denotes the derivative with respect to spectral parameter u. After a 
straightforward calculation, the equation ( p.lO|) does provide us with the expression 
( p. ID , whereas the second conserved quantity (|2.11|) has the form 



Here we would like to stress that both the Hamiltonian (2J.) and the conserved quantity 
( p.ll| ) should be understood as a global operators. It is meant that [H, J] = rather 
than Jj(j+i)(j_|_2)] = 0. The mutual commutativity of H and J convinces us of 

the existence of a quantum i?-matrix associated with the model (p.l|) . We shall present 
a rigorous proof of the integrability of the model in the next section. 

3. Integrability of the model 

It has long been clarified that the existence of the quantum i?-matrix which fulfills 
the Yang-Baxter relation is desirable for constructing integrable quantum chains. This 
suggests to us a way to verify the integrability of the model presented above. Indeed, 
following the paper [|ll|, we, after a cumbersome algebraic calculation, can find a class 



of solution to the Yang-Baxter relation 



V 



R (m, v)Loj{u) (g) Loj{v) = Lojiy) ® Loj{u) R {u, v), (3.1) 
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which is given as 



R (u,v) 
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with the Boltzmann weights 

Pi — {cosu cosv e' + sinT; sinu e~')p2, 



P4 = 
P9 = 

Pw = 
Pt- 

P3 = 
P5 = 



Pt 



{cosu cosv e~' + sinv sinu e')p2, 
(sin 14 cosve~'' — sinv cosu e')p2, 
(sinrt cosv e' — sinv cosu e~')p2, 



e p2, P2 =e P2, 

{cosu cosv e' — sinv sinu e~') 



cos^ u — sin^ V 



(cos It cosv e ' — sinv sin it e'] 



cos^ u 



sin^ V 



P2, 



P2, 



(cos u sin we' — sin v cos v e') 



cos^ u 



sin^ V 



P2, 



Pe 



cos u sm -u e — sm v 



cosv e 



cos^ u — sin^ V 



-P2 



and 



P8 = P3 - Pi, 
P7 = P5 - P4, 

I — h{u) — h{v), i 
which enjoy the following identities: 

P4Pi + PqPio = 1, 

PlP5 + P3P4 = 2, 

Pe Pe = P3P5 - 1- 



smh2h{u) _ U 
sin 2u 2 



(3.2) 



(3.3) 
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This i?-matrix with more distinct Boltzmann weights is indeed different to the one for 
the standard Hubbard model [1^, |Tl| and a twisted version which is associated 
with the Hubbard model with chemical potential terms. Running a Maple program we 
may check that the i?-matrix satisfies the Yang-Baxter equation 

Ruiu, v)Ri3{u, w)R23{v, w) = R23{v, w)Rn{u, w)Ri2{u, v). (3.4) 

So far we have built up the QISM machinism for the alternative Hubbard model and 
concluded the integrability of the model as well. On the other hand, a fermionic 
model is always favourable in the study of the condensed matter physics due to the 
clear distinguishment between the fermionic degrees of freedom and bosonic degrees of 
freedom. By performing the Jordan- Wigner transformations |11, 21 1, one may obtain 
the Hamiltonian of a fermionic model which is equivalent to the Hubbard model (|2.1|) : 



7V-1 ^11 

H = -Y1 J2^4j+i)s^i' + + ^ X](%T - o)('^(i+i)i - o) (3-5) 



Above aj^ and ajs are creation and annihilation operators with spins (s =t or |) at site 
j satisfying the anti-commutation relations 



33 ss y 



(3.6) 
(3.7) 



and Ujs = aj 

requires that the graded Lax operator related to the Hamiltonian ( p.5|) 



aj^ajs is the density operator. The integrability of the fermionic model ( |3.5| ) 



C 



u 



\ 



ie ^("^flj-T^ii 



ie~''^^^aj^a,i^ 
ie'^(")a,t^,-| 
e-'^''^93ian 



e 

-ie '^'"^'aj^aji e ■''"'a-^gji le 
must generate the graded Yang-Baxter relation 

V V 

TZ {u,v)Coj{u) (g) Coj{v) = Coj{v) Coj{u) TZ iu,v) 
with the graded i?-matrix which is given by 

W R {u,v)W-^ 



(3.1 



9j]93l 



V 

7^ iu,v) 



where 



W 



a 



(\ 

-i 



\0 








and 



f3 



3S 



smw 



smw 






1/ 



icos-uni 



9js 



COS u — (cos u + i sin M)n 



(3.9) 
(3.10) 

(3.11) 
(3.12) 



with the grading P(l) 
by 



P(4) = 0, P(2) = P(3) 



1. The monodromy matrix is defined 

(3.13) 
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such that the transfer matrices 

t{u) = stroT(M) (3.14) 
commutes each other for different values of the parameter u. It can be verified that an 



expansion of the logarithm of the transfer matrix ( |3.14| ) in powers of u will lead to the 



Hamiltonian (pj.5D as well as higher conserved quantities. 



We would like to remark that the model possesses the 5*0(4) symmetry if we 
consider a new representation of the r^-pair SU (2) algebra 

1=1 i=l 

and the C-pair spin SU (2) algebra 

L ^ L 

i=l i=l 

which comprise the SO {4) algebra. Taking into account the globality of these operators, 
one may show that the Hamiltonian (|3.5| ) commutes with the generators of the above 
two SU (2) algebras. This symmetry could be expected to complete all eigenstates of 
the Hubbard model like the case in the standard Hubbard model. Here the r^-pairing 
raising operator creating a pair of electrons with opposite spin on different sites could 
be interpreted as a delocalized Cooper pair. 

4. Algebraic Bethe ansatz 

Towards an exact solution of an integrable model, the algebraic Bethe ansatz seems to 
have more utility than the coordinate Bethe ansatz because the former not only provides 
us with the spectrum of all conserved quantities, but makes a close connection to the 
finite temperature properties of the model. There have been a lot of papers devoted to 
the study of the nested algebraic Bethe ansatz |2^ for the multistate integrable models 
with Lie algebra ( or Lie superalgebra ) symmetry. Following the so called ABCDF 
approach to solve the Hubbard-like models [|r^, we shall formulate the algebraic 
Bethe ansatz for the model in that which follows. To this end, as usual, we have 
to perform the ansatz step by step. However it is not necessary to restate all of the 
calculations used in solving our model because of the similarity to the routine proposed 



in the paper [|13 



In order to carry out the algebraic Bethe ansatz for this Hubbard model, we first 
need to find the eigenvalues and eigenvectors of the transfer matrix (|3.14| ) 



r|<l>„) = A|<l>„) (4.1) 
Following the prescription in [|13|, the eigenvectors of the transfer matrix are given by 
|$„) = $„.^|0). (4.2) 
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where the components of JF are coefficients of an arbitrary hnear combination of vectors 
$n and |0) is the pseudovacuum state, chosen here as the standard ferromagnetic one 

|0) = ®f=i|0)„ (4.3) 

where 



|0). 



1 




1 





(4.4) 



which corresponds to the doubly occupied state. We write the monodromy matrix T{u) 
( KWj ) as 

f B{u) Biiu) Biiu) F{u) \ 

Cr{u) A,^{u) Ai2(m) Er{u) 

C2{U) A2l{u) A22{U) E2{U) 

\ C^{u) C^{u) C^{u) D{u) J 

such that the necessary commutation relations between the diagonal fields and the 
creation fields can be derived from the Yang-Baxter algebra 



Tin) 



(4.5) 



T^i2{u,v) T (n) T (v) =T (v) T {u)nu{u,v). 



(4.6) 



above 

nu{u,v) = VR{u,v). 

Here V is the graded permutation operator. Let us ffist display an important 
commutation role, which reveals us a hidden nesting structure and the symmetry of 
eigenvectors, 

i 



B{u)B{v) = ^^^^^B{v)B{u).r{u,v) + 



piiu,v) 



+ 



p8{u,v)pi{u,v) 



F{v)B{u)^,{u,v) 



-F{u)B(v)Uu,v) 



(4.7) 



where 



r(u,vi 



with 



^iiu,v)= (0, /i(u,w), /2(M,t;),0) ; (2iu,v) = {O, p^{u,v), Pq{u,v),0) 
/I 0\ 

a{u,v) b{u,v) 
c{u,v) d{u,v) 

Vo 1/ 

fi{u,v) = Pq{u,v)p8{u,v) - p^{u,v)p5{u,v) 
f2{u,v) = Pq{u,v)p5{u,v) - p+{u,v)ps{u,v), 
P3{u,v)psiu,v) - p^{u,vf 



a{u, v) 
d{u, v) 



Pi{u,v)ps{u,v) 
ps{u,v)psiu,v) - peiu,vf 



b{u, v) = c{u, v) 



Pi{u,v)ps{u,v) 

ptiu,v)pQ{u,v) - p8iu,v)pr{u,v) 

pi{u,v)p8{u,v) 
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It turns out that the auxihary matrix f{u, v) is nothing but a gauged rational i?-matrix 
of an isotropic six-vertex model. If we adopt the parameterization introduced in [|T^ or 
2^ , explicitly, 



^ ^ _ smx ^ cos^^^Mx)^ ^^^^^^ (4.9) 
cos X sm X 



one may find that 



a(M,t')= — -: z d{u,v 



b{u, v) = c{u, v) 



u — v — U u — v — U' 

u — V 



u — v — U 



with 



-e{u,v) _ cosvsinu 



sm V cos u 

We shall see that the f-matrix ([4.8|) is related to the one of the isotropic six-vertex 
model via a proper gauge transformation, which does not change the spectrum of the 
spin sector. This seems to provide a new version of the i?-matrix, which does not have 
the difference property, for the isotropic six-vertex model. In view of the commutation 
relation (|4.7| ), the creation operators Ba-, Ea do not interwine. So it is reasonable that the 
eigenvectors of the transfer matrices are generated only by the creation operators Ba{u) 
and F{u). Following the argument in the paper [^, we may find that the n-particle 
vector can be determined recursively by the following relation 

$n(f 1, ■ ■ ■ , fn) = B{vi) ® $„-l(f2, ■ ■ ■ , fn) 

^ ip8(^^l,fi) l^. ip9{Vk,Vj) L 

F{Vi)^n-2{V2, ■ ■ ■ , Vj^i, Vj+i, ■ ■ ■ , Vn)B{Vj)\ --fk,k+l{Vk, ^^#.10) 

k=2 P^^'^^i "^i' 

Explicitly, the two-particle eigenvector reads 

$2^, V2) = B{v^) ® B{V2) + V2) ® F{vi)B{v2)—^ r. (4.11) 

^From the commutation relation ( |4.7| ), we can conclude that ■ ■ ■ satisfies an 

exchange symmetry relation 

^n{Vl, ■ ■ ■,Vj,Vj+i, ■■■,Vn) = j^-^ ' ^^"^^ | $n K , " " " " ■ , Vn) -f jj+i{Vj , Vj+i){4:.12) 

Pl[Vj, Vj+i) 

based on the following identity: 

Pi{Vj,Vj+i) 



, X / X / r6(^i+i>^i)-^(^i>^i+i) = 7 T6(^^j,^^i+l)-(4.13) 

Pl(^i+l>^i)P8(^^i+l,fi)Pl(^^i,^^i+l) p^[Vj,Vj^r) 

In above expressions, ^ plays the role of forbidding two spin up or two spin down 
electrons at same site. Also, F{u) creates a local hole pair with opposite spins. In order 
to manipulate the eigenvalue of the the transfer matrix ( p.l4|) we need the commutation 
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roles involving the diagonal fields over the creation fields. After some algebra, from the 
Yang-Baxter relation ([4 .61 ) we have 

1 



Biu)Baiv) = Pl^:^Ba{v)B{u) — 

D{u)B^{v) = 'P^^p^B^{v)D{u) - — ^ 

P8{U,V) P8{V,U 



-Baiu)B{v).rii{v,u), 



(4.14) 



P8{U,V) 



Uu,v). {E*{u)(g)A{v) 



(4.15) 



Aat{u)Ba{v) = ^^Biv)^Aiu).fiu,v) 



p9iu,v) 



p9iu,V 



-fj2{u, v).B{u) ® A{v) 



+ 



pg{u,v)p8{u,v) 

+ / , E*{u)B{v) ® i2{u, v). 



(4.16) 



Above we introduced the notations 

' pt{u,v) 



Aiu) 
B 

& 



P2(«>^) 

p^(u,t;) 

pt{u,v) 

An(u) Ai2(m) 

^21 (m) ^22(m) 



(i?l,52), C 



(C4, C5) , E* 



C^2 

^1 
E2 



(4.17) 



In order to determine the eigenvalue of the transfer matrix ( |3.14| ) acting on the mult- 
particle eigenstates we need to consider the commutation relations for the creation field 
F{u\. 



B{u)F{v) 



^^("'")-F(.)5(«) + ^4^F(«)i?(.) 



D[u)F[v) 



Ps{v,u) 
i 

psiv^u 
i 



v,u} 



B{u)®B{v) .iliv^u), 



(4.18) 



F{v)D{u) + ^^^j^F{u)D{v) 
Ps[v,u) 

iiiy.u). \b{u)®b{v) 



A{u)F{v) 



ps{v,u 

p9{u,v)pio{u,v) 



(4.19) 



F{v)A{u) 



pg{u,v)pio{u,v) 



fi2{u, v).F{u)A{v).fi2{u, v) + -. — 7'n2{u, v).B{u) (g) E*{v 

lp9{u,V) 
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^ -E*{u)^B{v).fi2{u,v), (4.20) 



ipio{u,v 



B{u)F{v) = 'P^p^F{v)B{u) + — ^r/2(«, v).B{v)F{^^), (4-21) 
pi{u,v) pi{u,v) 

F{u)B{v) = -'-Pl^f^B{v)F{u) + — L^r}i(M, v).Fiv)Biu). (4.22) 
pi{u,v) Pi{u,v) 

Finally, if we adopt the variables z±{vi) used in ,i.e. 

sm Vi cos Vi 

and make a shift on the spin rapidity = Xj + U/2, the eigenvalue of the transfer 
matrix ( p.l4| ) is given as (up on a common factor) 

sinM(l + z_{vi) / z+{u)) 



t{u) I ^^{vi,---,Vn)) = < [Z-iu)]''Yl 



^J- COS'u(l — Z_{Vi) / Z-{u)) 



-\-[Z^{U \ I I r ——- 

cos-u^l — Z_[Vi)Z^{u)) 



-TT sinM(l + z_{vi)/z+{u)) T-r {u ~ Xj + u/2) 
l\ cosm(1 - z.{v,)/z^{u)) l\ {u-Xi - u/2) 

A sin^(l + z.iv,)z.{u)) fr ju-Xi- 3U/2) 1 



cos Ml 

«=1 



provided that 



l-e.)l^^nf4^ (4.25) 

7=1 (^^i - A; - U/2) 

A (A, - f). + [7/2) _ -|^^ (A, - A, + [/) 

(A, - - t//2) J-1 (A^._A,-f/) 

t i, 



where 



J = 1, ■ ■ - M, z = 1, ■ ■ ■ ,n. 

If we express the variable Z-{ui) in terms of the (hole) momenta ki by Z-{ui) = e^^^, 
from the relation (|4.9| ), the energy is given by 



E„ = -(Ar/2-ri)[/-^2cosfci. (4.27) 

■t=i 

Using the momenta ki instead of the charge rapidity Vi via the relation ( [4.9[ ) and making 
a scaling on the spin rapidity Xj as A-, = — | Aj, then the Bethe equations ( [4. 251 ) and 
i W2^ ) read 
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" (sinA;.-A,-f) _ ^ {X, - A, + f ) 

ii(sinA:.-A, + f) [\ (A,-A.-f)' ^^''^^ 

L i, 



j = !,■■■ M, i = l,---,n. 
5. Conclusions and discussion 

We have proposed an integrable variant of the Hubbard model with off-site Coulomb 
interaction. The integrability of the model was verified by showing that the quantum R- 
matrix satisfies the Yang-Baxter equation. It was argued that the model possess SO (4) 
symmetry, however, it contains a new representation of 77-pairing SU{2) algebra and 
C-pair spin SU{2) algebra. By means of the nested Bethe ansatz, we have presented 
the spectrum of the Hamiltonian, eigenvectors and the Bethe ansatz equations for the 
model with periodic boundary conditions. We found that the model exhibits a gauged 
r-matrix of the isotropic XXX model, which plays a crucial role in solving the model. 
Under periodic boundary conditions the alternative model and the standard Hubbard 
model share the same spectrum and Bethe ansatz equations. This result comes as 
no surprise since the model proposed here is obtained from the usual Hubbard model 
through the transformation 

Q| Ci|, di C(i+i)| (5.1) 

which also maps the standard SO (A) symmetry generators to ( p. 151 , ^1161) . However, the 



new quantum i?-matrix will lead to different open boundary conditions from the ones for 



the standard Hubbard model p5|, since the transformation (|5.1|) will not preserve the 
open chain. In turn, the differences in spectrum for the two models would be apparent 
in the case of open boundary conditons. This seems to open an opportunity to identify 



new boundary impurity effects |27, Bq, G9i in a Luttinger liquid. An interesting problem 



is to identify the boost operator for the spectral parameter extension of this new model. 



which can iteratively generate all of the conserved currents, using the results of |T^. We 
we shall be focusing on these problems in near future. 
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